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Abstract 

We calculate the shift in the critical temperature of Bose-Einstein conden- 
sation for a dilute Bose-Fermi mixture confined by a harmonic potential to 
lowest order in both the Bose-Bose and Bose-Fermi coupling constants. The 
relative importance of the effect on the critical temperature of the boson-boson 
and boson-fermion interactions is investigated as a function of the parameters 
of the mixture. The possible relevance of the shift of the transition temperature 
in current experiments on trapped Bose-Fermi mixtures is discussed. 



I. INTRODUCTION 

The achievement of Bose-Einstein condensation in ultracold, trapped dilute alkah gases, 
beyond reahzing a striking and spectacular experimental confirmation of a long-standing, 
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fundamental prediction of quantum theory, has initiated and stimulated as well a whole new 
field of research in the physics of quantum gases in confined geometries [1,2]. Nowaday, 
many different atomic species and isotopes can be succesfuUy cooled and trapped in the 
gaseous state, to investigate in exceptional conditions of purity and isolation the properties 
of interacting Bose and Fermi systems, or Bose- Fermi mixtures. 

In particular, the experimental realizations of trapped gaseous mixtures of bosons and 
fermions are both an interesting new instance of a quantum many-body system and a very 
useful tool to reach the regime of quantum degeneracy for a Fermi gas via sympathetic 
cooling of the fermions by the bosons [3-7]. From a theoretical point of view, dilute Bose- 
Fermi mixtures have been the object of recent investigations addressing, for example, the 
determination of the density profiles of the two components in trapped systems [8], the 
problem of stability and phase separation [9], and the effect of boson-fermion interactions 
on the dynamics [10] and on the ground-state properties [11] of the mixture. 

Boson-fermion interactions in a Bose-Fermi mixture can induce a net attractive interac- 
tion between the fermions, thus introducing a further mechanism toward the achievement 
of the BCS transition in trapped Fermi gases [12]. In this paper we address the reverse 
problem, i.e. how the transition temperature of Bose- Einstein condensation is affected by 
the presence of the fermions in a trapped mixture. The shift of the transition temperature Tc 
due to interactions in a pure trapped Bose system has been calculated within the mean-field 
approximation in Ref. [13]. In the present work we extend the perturbative methods of Ref. 
[13] to obtain the shift of Tc to lowest order in both the Bose-Bose and Bose-Fermi coupling 
constants. To this order the effect on Tc of boson-boson and boson-fermion interactions are 
independent and add linearly. The relative importance of the two effects depends on the 
relevant parameters of the trapped mixture: number of bosons and fermions in the trap, 
ratio of the masses and of the oscillator frequencies for the two species and the ratio of the 
Bose-Bose and Bose-Fermi coupling constants. The calculation is carried out in local density 
approximation which is valid provided that the number of bosons and fermions in the trap 
is large. Finite size effects have not been included in the present treatment. 
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The plan of the paper is as follows: in Section II we generahze the scheme derived 
for pure Bose systems in Ref. [13], to include the effects of boson-fermion interactions. In 
Section III we derive analytical results for the shift of Tc in the limits of a highly degenerate 
(Thomas-Fermi) and a classical (Boltzmann) Fermi gas, and we provide the full numerical 
solution for the intermediate regimes. We finally compare the theoretical predictions with 
the current experimental situations and we draw our conclusions. 



II. THEORY 

In a non- interacting Bose gas confined by the external harmonic potential Vgf^(r) = 
mBifjjIx^ + ijJyU'^ +uj1z^)/2, the critical temperature for Bose- Einstein condensation (BEC) 
is given by 

f N \^^^ 

keT^ = hujB i-^A OMhujB N'J' , (1) 

where cub = ('^I'^y'^z)^^^ is the geometric mean of the oscillator frequencies and tub, Nb are 
respectively the particle mass and the number of bosons in the trap. The above result is 
obtained using the local density approximation (LDA), where the temperature of the gas is 
assumed to be much larger than the spacing between single particle levels: fc^T ^ fiuJx,y,z- 
In this case the density of thermal atoms can be written as 

<(r) = (Af )-^^3/2(exp{-[Kf,(r) - A^bI/^bT}) , (2) 

where A-^ = h{2n / mskBTy^'^ is the boson thermal wavelength, and (73/2(2;) = J2^=i ^'^ / it-^^"^ 
is the standard Bose function of order 3/2. At T — the boson chemical potential takes 
the critical value hb — IJ'c — 0) corresponding to the bottom of the external potential, and 
the density n^(0) in the center of the trap satisfies the critical condition n^(0)(A^o)^ = 
C(3/2) ~ 2.61 holding for a homogeneous system. 

Finite size effects modify the prediction of the critical temperature (1) resulting in a 
reduction of T°. The first correction due to the finite number of atoms in the trap is given 
by [14]: 
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where Cob = (i^x + cjj; + ujz)/'i is the arithmetic mean of the oscillator frequencies. 

Interparticle interactions have an effect on the BEC transition temperature as well. The 
presence of repulsive interactions has the effect of expanding the atomic cloud, with a con- 
sequent decrease of the density. Lowering the peak density has then the effect of lowering 
the critical temperature. On the contrary, attractive interactions produce an increase of the 
density and thus an increase of Tc- This effect, which is absent in the case of a uniform gas 
where the density is kept fixed, can be easily estimated within mean-field theory. For pure 
bosonic systems the shift STc = T^. — T° has been calculated in Ref. [13], 

(§) = -1.33 <• . (4) 

to lowest order in the coupling constant Qbb = 4:Tih^aBB/'iTiB- In the above equation qbb 
is the boson-boson s-wave scattering length and Ib — ^Jh/fnBOJB is the harmonic oscillator 
length. Result (4) has been obtained within LDA and neglects finite size effects. 

In the case of trapped Bose- Fermi mixtures the shift of Tc, due to both Bose-Bose and 
Bose-Fermi couplings can be calculated in mean-field approximation using the methods of 
Ref. [13] . The transition temperature Tc of a trapped Bose gas is defined by the normalization 
condition 

Nb = j drnB{r,Tc,i^c) , (5) 

where is the thermal density of bosons and /^c is the critical value of the boson chemical 
potential. Within LDA the boson density above Tc is given by 

riBir) = (Af )-^^3/2(exp{-[K?/(r) - piB]/kBT}) , (6) 

where 

K//(r) = Kxt(r) + 2gBBnB{r) + gBFUpir) , (7) 
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is the effective potential acting on tlie bosons wfiicfi is generated by tfie external field V^^ and 
by the mean field produced by interactions with the other bosons and with the fermions. 
Notice the factor 2 present in the Bose-Bose contribution and absent in the Bose-Fermi 
term due to exchange effects. In the above equation npir) is the fermion density and 
qbf — I'Kfi^aBF /"^R is the Bose-Fermi coupling constant, fixed by the boson-fermion s-wave 
scattering length aBP and by the reduced mass ruR = mBrriF/ {ms + m^), where mi? is the 
fermion mass. 

For a fixed value of the boson chemical potential jiB and a fixed temperature T, the 
boson density (6) can be expanded to first order in qbb and qbf as 

nB(r, T, fiB) = n%{r, T, i^b) - [2^sB<(r) + ^^^^(r)]-^ , (8) 

djiB 

in terms of the non-interacting boson (2) and fermion density 

<(r) = (A|;)-V3/2(exp{-[K^,(r) - /.^j/fc^T}) . (9) 

In the above equation A-^ = h{27r / mpkBTy^'^ is the fermion thermal wavelength, and f3/2{x) 
is the Fermi function of order 3/2 defined as 

/3/2(x) = ^/ dz^^. (10) 

Result (9) has been obtained in LDA for a Fermi gas in the trapping potential V^^{v) — 
mF{u'^x^ ^-u'^y^ -\-uj'^ z^) /2. The fermion chemical potential /xp is fixed by the normalization 
condition 

Nf = j dr n%{r) , (11) 

where Nf is the total number of fermions in the trap. The condition of validity for LDA 
requires the Fermi temperature of the fermionic system to be much larger than the harmonic 
oscillator energies /cb^f ^ ^^x,y,z- ^ non- interacting trapped Fermi system the Fermi 
temperature, or equivalently the Fermi energy, is given by /c^Tp = ep = hup^GNFy^^ , where 
ujf — i^'x^'y^zy^^ is the geometric mean of the fermion oscillator frequencies. 
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To first order in qbb and qbf, the critical value //c of the boson chemical potential can 
be written as 

Ate = + 2gBBnUr = 0) + ^BF<(r = 0) . (12) 

By writing — T° + 5Tc, one can expand Eq. (5) obtaining the following result for the 
total relative shift of the condensation temperature: 

^^(^JIA +(^] '^SBB I dr dn%/diiB[n%{v = 0) - n%{v)] 

7? WJbb \ T'c J BF TO J dv dn%/dT 

Qbf I dr dn%/dnB[n%{r = 0) - n%{r)] 
TO Jdrdn%/dT ' ^ ^ 

where the derivatives of the non-interacting boson and fermion densities and n^, are 
evaluated at the ideal critical point /i^ — 0, T — T°. The first term {STc/T^)^^ in the above 
equation accounts for interaction effects among the bosons and coincides with the shift (4). 
The second term {6Tc/T^)^p accounts instead for interaction effects between bosons and 
fermions, and its determination will constitute the main result of the present paper. Some 
comments are in order here, (i) The shift ST^. derived above is a mean-field effect which 
originates from the fact that in a trapped Bosc-Fcrmi mixture the total number of bosons 
and the total number of fermions are fixed, but not the density profiles of the two species. 
This effect is peculiar of trapped systems, since it vanishes identically in the case of uniform 
systems, and should not be confused with the shift of occurring in homogeneous Bose 
systems, which is instead due to many-body effects [15]. (ii) The shift originating from 
the Bose-Fermi coupling, similarly to the one arising from the Bose-Bose one, is negative if 
Qbf > and is positive if qbf < 0. If obb and a^F have opposite sign, the corresponding 
shifts of Tc go in opposite directions, (iii) Result (13) holds to lowest order in qbb and qbf 
and, since it has been obtained using LDA, is exact if the number of bosons and fermions 
is large. Finite-size corrections are not included in (13). For a finite system, a reasonable 
estimate of the total shift of the critical temperature can be obtained by adding to result 
(13) the finite-size correction (3) of the non- interacting model. 
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III. RESULTS 



We now concentrate on the relative shift {5Tc/T^)^p due to the boson- fermion interac- 
tion. First of all we observe that 

^ - (A5)b? ^.Mexp[-Kf,(r)/;=BTa , (14) 

and r°/ dr dnP^/dT — SNb, where the derivatives are evaluated at the condensation point 
of the non-interacting gas /J^ — 0, T — T^. Using Eq. (9), the relative shift can then be 
rewritten as: 

/ STA _ _9bf_ 1 

X [/3/2(exp{/.^/fc5Te}) - /3/2(exp{[/.^ - K^,(r)]/fcsT°})] . (15) 

In the following we shall assume that even if the trapping potentials of bosons and 
fermions can have different oscillator frequencies, nevertheless LO^jio'^ = ujylLo'y — oJz/oj'^ — 
cub/^^f, i-e. the anisotropy is the same for the bosonic and fermionic trapping potentials. 
This is always the case in today's experiments, and assuming otherwise would introduce 
unnecessary complications. In fact, the assumption of equal anisotropics holds in general 
in magnetic traps since the confining potentials depend only on the (common) external 
magnetic field, the magnetic moments, and the masses of the atoms. Eq. (15) contains the 
fermion chemical potential ij,f{Nf, T°) which has to be determined from Eq. (11). Eqs. (15) 
and (11) have then to be solved simultaneously. We notice that Eq. (11) can be rewritten 
in dimensionless form as 



oo +2 

T| = 3 / dt ^ , (16) 

^ exp(t-/iF) + l 







where we have introduced the reduced chemical potential flp — I^'f/^bT^ and the reduced 
Fermi temperature Tp — Tp/T^. Eq. (16) reveals that /ip is only a function of Tp, which 
in turn is a measure of the degeneracy of the Fermi gas at T = T°. In terms of jlp and Tp 
Eq. (15) then becomes 



X j ds s"^ gi/2{exp{—TFas^}) 

X [/3/2(exp{/iF}) - /3/2(exp{/iF - Tf s^})] . (17) 

In writing Eq. (17) we have rescaled each integration coordinate by the appropriate Thomas- 
Fermi radius of the fermion cloud i?- = {2eF/mF0j'-^Y/'^. We have then introduced the mean 
Fermi radius Rf — {R'^R'yR'^Y^^ and named a — mBOj\lmFOJ^- Since jiF depends only on 
Tf through Eq. (16), the integral in Eq. (17) above depends only on the values of the two 
parameters Tf and a. 

The system of Eqs. (17) and (16) for general Tf and a can only be solved numerically, and 
later we shall present the full numerical results for some specific choices of the parameters. 
However, analytical solutions exist in two limits: when Tp ^ 1 (i-e. Tf ^ T^) where the 
Fermi gas is completely degenerate at T = T° (Thomas- Fermi regime), and when Tp ^ 1 
(i.e. Tf -C T^) so that at T° fermions behave as a classical gas (Boltzmann regime). 

In order to clarify the connection between the two limits and the general numerical 

solution, it is useful to further manipulate Eq. (17). By explicitly evaluating the prefactor, 

it can be finally recast in the convenient form 

(6TA 25/3 
\T^J^^ 3V6;rC(3) \mB 



f^^ + l) ^iVr •F(rp,a), (18) 



where 



F{fF,a)^a'/''fF J ds g'i/2(exp{-TF a s^}) 

X [/3/2(exp{/iF}) - /3/2(exp{/iF - Tf s^})] , (19) 



and £f = yh/rriFiOF is the fermionic oscillator length. Notice the formal analogy between 
Eq. (18) and Eq. (4) for the shift {8Tc/T^)bb due to the boson-boson interactions alone. 

Let us begin by considering the Thomas- Fermi limit (Tp ^1). In this limit the chemical 
potential of the fermions /xp tends to the Fermi energy ep = ksTp- Thus fiF — Tp ^ 1- 
The limit of the Fermi functions in Eq. (10) for x — > oo is fz/2{x) ~ 4(lnx)3/^/3y^. This 



implies that the density profile of the fermion cloud takes the well known Thomas- Fermi 
shape 

n%{v) = 4(0) [l - {x/R'^r - {y/R'^f - {z/R',ff' , (20) 

with n%{0) — (2ei?mF/^^)^/^/(67r^), whenever the expression inside the square brackets is 
positive, and n%{T) = otherwise. 

The function F{Tf, a) then goes to the limiting form 

X e-"^^°^'[l-(l-s2)3/2] , (21) 

since 5^1/2(2;) = E^=i 

We obtained Eq. (21) in the limit Tp ^ 1. Therefore, if a is not too small (so that 
Tpa ^ 1 still holds), then, for every n in the series, the exponential is non- vanishing 
only for values of s <S 1, and we can adopt the expansion 1 — (1 — s^)^/^ ~ 3s^/2. The 
integral in Eq. (21) becomes Jq dss^e""^^"*^ ~ 30r/[8(nTF a)^/^]. Finally, therefore, 
F{Tf, a) — > 3C(3)/4ci; and the Thomas- Fermi prediction for the relative shift reads 



where 3^/^/ (2^/^7r) ~ 0.304. We notice that in the Thomas-Fermi regime the shift is inde- 
pendent of the number of bosons Nb and varies as the first inverse power of the parameter 

a — tub^j'^b/^f^'f- 

We now consider the Boltzmann limit for the Fermi gas (Tp -C 1). In this case the 
chemical potential is large and negative and depends on Tp as: jxp ^ 1^{(^f)^/6}. In 
the limit x — > 0, h/2{x) ~ x, and then 



F{Tp,a)^ LZ^^ / dss' 

6 n^'-^ Jo 



X 



(23) 



Evaluation of the integrals is straightforward and yields 



F{f^,a)^^{fFf'-f{a), (24) 

with 

°° / 1 1 \ 

/(«) = E - - i/2( ^ -n3/2 ' (25) 

so that: 

where the numerical prefactor is ~ 0.025. 

We notice that f{a) in Eq. (25) is a monotonically decreasing function of a. In particular, 
one finds the following behaviours: f{a — > 0) = 7r^/6, /(I) ~ 0.85, and f{a — > oo) = 
3C(3)/2q;. As one should expect, in the Boltzmann limit (Tp <g 1), the Bose-Fermi shift is 
negligible. 

We now turn to the full numerical solution of Eqs. (18), (19), and (16) for more general 
values of the degeneracy parameter Tp. In Fig. 1 we show the dimensionless function 
F{Tf, a) as a function of Tp for three different values of the parameter a = mBU!%/'mFUJF, 
a — 0.1, 1, and 10. For fixed a, F{Tf, a) is a monotonically nondecreasing function of Tp, 
which saturates for Tp ^ oo at the value predicted in the Thomas- Fermi regime 3C(3)/4q; ~ 
0.9 a~^. For fixed Tp, F{Tf, a) increases by decreasing a. For the largest value of o; (a = 10) 
the function F reaches its asymptotic Thomas- Fermi value already at Tp ~ 5. For a — 1 and 
a = 0.1 the function saturates for larger values of Tp not shown in the figure. The reason 
for this difference can be understood by recalling that the Thomas-Fermi result requires not 
only Tp » 1, but also Tp » (see the discussion below Eq. (21)). 
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FIG. 1. Dimensionless function F{TF,a) as a function of Tp for the values a = 0.1 (dotted 
line), a = 1.0 (dashed line), and a = 10 (solid line). 

The physically relevant regimes in current experiments fall roughly around a ~ 1 and 
Tp — 1- In this respect, a particularly interesting situation is the one realized in the Florence 
experiment [7], where a quantum degenerate trapped atomic mixture of fermionic ^''K and 
bosonic ^^Rb has been recently produced. One of the appealing features of this system is 
that the measured boson-fermion scattering length is large and negative: qbf — —22 nm, 
giving rise to a fairly strong attractive boson-fermion interaction. The shift {6Tc/Tj^)^p is 
thus positive and opposite to the shift {STc/T^)^^, since for pure ^^Rb the boson-boson 
scattering length is qbb — 6 nm, giving rise to a repulsive boson-boson interaction. In the 
Florence experiment the two atomic species are magnetically trapped, and are both prepared 
in their doubly polarized spin state. These states experience the same trapping potential so 
that a — ■mBUJ%/'mFUJF — 1, while the number of bosons and of fermions are respectively 
A^B = 2 X 10^ Nf = 10^ so that Nf/Nb = 0.5, and Tp = Tp/T° ~ 2.3. For the conditions 
of the Florence experiment the shift (4) due to the boson-boson coupling turns out to be: 
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{5Tc/T^)qq ~ —0.037, and is comparable with the shift (3) due to finite size effects, which 
is given by: {6TJT^)^^ = -0.044. For a = 1 at 7> ~ 2.3 the function F is at about 
1/3 of its asymptotic value in the Thomas- Fermi regime, resulting in a Bose- Fermi shift 
considerably smaller than the Bose-Bose one: {5Tc/T^)^p ~ 0.012. In Fig. 2 we show the 
shift {6Tc/T^)^p as a function of the ratio Nf/Nb, with all the other parameters entering 
Eq. (18) fixed at the values of the Florence experiment [7]. In the same figure we include 
as a reference value the modulus of the boson-boson relative shift | {STc/T^)bb h calculated 
using the values of the parameters given by the Florence experiment. 

0.045 I 1 1 1 1 1 1 1 1 1 1 




1 23456789 10 



Nf/Nb 

FIG. 2. Boson-fermion relative shift {5Tc/T^)bf from Eq. (18) (solid line) as a function of the 
ratio Nf/Nb- Horizontal dashed line: value of the modulus | (STc/T^) \ of the boson-boson shift 
(4). All other parameters, except the number of fermions Nf, have been fixed at the values of the 
Florence experiment. 

From Fig. 2 we see that, while in the present experimental situation the boson-fermion 
shift is about 1/3 of the boson-boson one, by increasing the number of trapped fermions the 
two shifts become comparable at Nf — 5Nb- The boson-fermion shift is instead dominant 
at still larger values of A^^. It is important to remark that, even if the Bose-Fermi shift of 
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the critical temperature is a small effect for the present experimental conditions, it might 
be observable. Since the fermions can be eliminated from the trap, one can look for the 
differences in the transition temperature with and without fermions. 

In conclusion, we have determined the relative shift of the critical temperature of Bose- 
Einstein condensation in a trapped atomic Bose-Fermi mixture to lowest order in both the 
boson-boson and the boson-fermion coupling constants. We have determined numerically the 
general behaviour of the boson-fermion shift, and we have provided full analytical solutions 
in the quantum degenerate Thomas- Fermi regime and in the classical Boltzmann regime. We 
have applied our predictions to a specific experiment (the Florence experiment [7], chosen 
for the interesting value of the Bose-Fermi scattering length), and discussed the relative 
importance of the shifts due to boson-boson and boson-fermion interactions. 
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